DTP/93/12 
RAL-93-014 
March 1993 



Phenomenological constraints on the flavour 
asymmetry of the nucleon sea 



A.D. Martin and W.J. Stirling 



\Q | Department of Physics, University of Durham 

Durham DH1 3LE, England 

> \ and 

^ ! R.G. Roberts 

^ ' Rutherford Appleton Laboratory, 

m ■ Chilton, Didcot 0X11 OQX, England 

a,: 

£ni Abstract 
r- ! 

We study the possible flavour asymmetry, u ^ d, of the light quark sea 
distributions of the proton. We discuss the information that is at present 
available from data on deep-inelastic lepton-nucleon scattering and from Drell- 
Yan production on various nuclear targets. We show that the ratio of dilepton 
yields on hydrogen and deuterium targets is very sensitive to u — d. 



The increased precision of deep inelastic lepton-nucleon scattering data has con- 
siderably improved our knowledge of the parton distributions of the proton. Until 
recently all global structure function analyses assumed that the light quark sea dis- 
tributions satisfied the equality u — d. However, motivated by the improvement of 
the data, particularly at small x, and by the consequences for the Gottfried sum rule, 
analyses [|T], |2|, [3] have been performed with u ^ d. Here we wish to explore how large 
an inequality the available data can tolerate and to see how well future experiments 
can pin down u — d. 

Valuable insight into the problem is obtained by expressing (to leading order) the 
observed muon and neutrino deep inelastic structure functions in terms of the parton 
densities 



where N is an isoscalar nuclear target, and where we have included the s quark, but 
neglected the (small) c quark, contribution. These four observables determine four 
combinations of parton densities, which we can take to be u + u, d + d, u + d and s. 
In other words data at a given x, Q 2 cannot verify whether u = d or not. 

The above discussion is very simplistic because in practice we must (i) include 
the next-to-leading 0(a s ) contributions, which require knowledge of the gluon dis- 
tribution, (ii) allow for possible adjustments in the overall relative normalisation of 
the data sets, (iii) include the heavy target corrections of the neutrino data, and 
(iv) include screening corrections to allow F% n to be obtained from deuterium data. 
Assuming that this has all been done correctly, we compare in Fig. 1 the parton 
distributions at Q 2 = 20 GeV 2 from the two most recent global analyses Ejfl. If 
we make a detailed comparison at, say, x = 0.15 then we see significant differences 
between u, d, u, d from the two sets of partons. For CTEQ || we have d ~ 2u, while 
for MRS Q we have d ~ However this difference is compensated by w(CTEQ) 
> w(MRS) and d(CTEQ) < d(MRS) so that the observable combinations u + u, d + d 
and u + d remain essentially the same in both analyses, as indeed they should be for 

lr This difference is a consequence of the flexible (and independent) parametrizations of the u, d, s 
distributions allowed in the CTEQ analysis || . However the most dramatic difference between the 
CTEQ and MRS distributions is in the behaviour of the s quark distribution (the dash-dot curves 
on Fig. 1). The freely parametrized CTEQ s distribution becomes relatively large at small x to 
accommodate the F£ N data for x < 0.1. On the other hand MRS [Q, || adopted a less flexible 
parametrization, with s = j(w + d) at Q 2 — 4 GeV 2 , and allowed for uncertainties in the heavy 
target corrections to the neutrino data; the (factor of 2) suppression of the s quark distribution was 
chosen so as to be in better accord with dimuon production in deep inelastic neutrino scattering [H . 
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values of x where accurate data exist. So, as mentioned above, without additional 
input, the structure function data on their own do not test u ^ d. 
The motivation for u 7^ d comes from the NMC measurement [[J 

r - 8 dx 

/ — (F 2 V — F 2 n ) = 0.227 ±0.007(stat.)±0.014(sys.) (5) 

J0.004 X 

at Q 2 = 4 GeV 2 , as compared to the Gottfried sum rule 

Igsr= / — (Ff-FD = ~ dx{u v -d v ) + - dx{u-d) 
Jo x 6 Jo o Jo 

= - if u = d , (6) 

where here u and d have been expressed as the sum of valence and sea distributions, 
e.g. u = u v + u s and u = u s . A straightforward comparison of (5) and (6) implies 
that d > u, and indeed from the lack of Regge f 2 — a 2 exchange degeneracy we expect 

d-u = Ax- aR (l-xy s (7) 

with the small x behaviour governed by the Regge (meson trajectory) intercept ~ 
0.5. In the global analyses of refs. ]T|, we parametrized d — u in this way and found 
Igsr — 0.24 — 0.26. However it is possible to achieve an equally acceptable global 
fit of all the structure function data with u = d (and hence Igsr = |), albeit with a 
contrived small x behaviour of the valence distributions ffl]. 

Fig. 2 shows the ingredients of the Gottfried sum rule for various sets of parton 
distributions with u 7^ d. The central part of the figure, which shows F 2 P — F 2 n on 
a logarithmic x scale, is a visual display of the integrand of the sum rule at Q 2 = 7 
GeV 2 . The data are from NMC 0. The top part of the figure shows the accumulated 
contribution to the sum rule as we integrate down from x — 1; thus the limiting values 
of the curves as x — > give Igsr of (6). The central curves in the shaded bands are 
obtained from the MRS(Dq) partons with d — u parametrized at Q 2 = 4 GeV 2 as in 
(7); these partons give Igsr = 0.256 ||. 

In order to explore the allowed variations in the shape of the x behaviour of d — u 
we have repeated the global analysis of ref. || with the difference having the extended 
parametrization 

d-u = Ax~ aR {l -x) vs (l + 7x + 5x 2 ), (8) 

with A chosen to maintain the normalisation Igsr — 0.256. First we set 5 = 0. We 
find that we can maintain the quality of the fit to the data, that was achieved by the 
Dq fit of ref. 0, provided that 7 lies in the range —8 < 7 < 32. These extreme "hard" 
and "soft" versions of d — u of Dq correspond to the boundary curves of the shaded 
regions shown in Fig. 2. We may thus regard the shaded region as a first qualitative 
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estimate of the uncertainty in d — u. We see that this range of fits only modifies 
F% v — F% n for x < 0.05; indeed deep-inelastic data impose only a weak constraint on 
d — m at larger values of x. For a determination of d — u at larger x we should compare 
Drell-Yan production of pp and pn origin. A hint of what may be possible is shown 
in Fig. 3, which compares Drell-Yan data |7j obtained from a neutron rich (tungsten) 
target with those from isoscalar (deuterium and carbon) targets. We discuss such 
Drell-Yan comparisons, in detail, below. 

We have also repeated the global analysis of the deep inelastic data with 5^0 
in the parametrization (8) of d — u, in an attempt to reproduce the qualitative trend 
of the Drell-Yan data of Fig. 3. The predictions of a fit with 7 = —60, 5 = 300 
are shown by the dashed curves on Figs. 2 and 3. These have a radically different 
shape and, in fact, the quality of the fit to the deep inelastic data diminishes; we 
regard this as an example of extreme behaviour of d — u. For completeness we show 
by dot-dashed curves the predictions of the CTEQ1M analysis, which is in the spirit 
of the Dq fit but with u, d, s much more freely parametrized. We see a dramatically 
different behaviour at small x; indeed F% p — F% n becomes negative for x < 0.006 and 
as a result the value of Igsr is well below 0.2. 

The distributions shown in the lower part of Fig. 2 therefore offer (i) a set of 
conservative variations of d — u of the Dq fit which are consistent with all the deep 
inelastic data, together with (ii) two examples of unusual (but interesting) behaviour 
which arise if the sea distributions are more freely parametrized. 

The Drell-Yan process provides a complementary method of investigating the 
antiquark content of the nucleon. In leading order, for a proton beam on a nuclear 
target A, 

d 2 a Aira 2 1 ^-^ 2 r , s _. . . , ... 

dAPd^ = 9APS x~^x~ A ? C « fc^M + ^« ' ^ 

with t = M 2 /s = x p xa and xp = x p — xa- The crucial point is that a comparison 
of the cross sections on targets with different numbers of protons and neutrons leads 
directly to information on d—u. Two methods have been suggested. The first relies 
on the fact that at large xf the cross section (9) is dominated by the annihilation of 
u quarks in the proton with u quarks in the target. If the latter contains Z protons 
and A — Z neutrons, then 

In particular the ratio of non-isoscalar to isoscalar target cross sections is [0 

p fV> ~ 1 + fl 2Z \ d{x)-u{x) 
Ra/Is(x) ~ 1 + l 1 " Tj d(x)+u(x) ' 
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with x = xa- The ratio for tungsten (with 1 — 2Z/A = 0.195) to isoscalar targets has 
recently been measured by the E772 collaboration 0. A comparison of their data 
with the predictions of the various parton sets discussed above is shown in Fig. 3. 
The open circles at small x show the data before a correction for nuclear shadowing 
is applied |]]. It turns out that the xp values are largest for the data points at small 
x, and it is only for these that the approximation (|TT| ) is valid. At larger x, the 
smaller xp values lead to a contamination of the u(x p )u(xa) annihilation with the 
other processes, and a corresponding decrease of the sensitivity to d — u. The shaded 
band in Fig. 3 again corresponds to the range 5 = 0, —8 < 7 < 32 of the d — u 
parametrization of (8), evaluated using the full Drell-Yan next-to-leading order cross 
section. The Dq prediction (solid line) is roughly in the middle of this band. The 
Sg prediction (not shown) has R ~ 1 for small x, falling slightly to R ~ 0.99 at the 
larger x values. All these predictions are clearly consistent with the measured ratio. 
Notice that there is perhaps some evidence that R < 1 (i.e. d < u) for x < 0.15. As 
mentioned above, it is not difficult to derive a d — u form which follows this trend, as 
shown by the dashed line in Fig. 3 which corresponds to 5 = 300, 7 = —60, though 
this parametrization leads to a poorer overall fit to the deep inelastic data. Note also 
that the CTEQ1M prediction (dot-dashed curve), which has d — u large and positive 
for the range of x spanned by the data, appears to be slightly disfavoured. 

A much more powerful Drell-Yan measurement is to compare the cross sections 
on proton and neutron (in practice hydrogen and deuterium) targets ||. Two such 
experiments have recently been proposed || [H|. Fig. 4 shows the next-to-leading 
order QCD prediction for the asymmetry 

Ady = <?(PP) - <?(pn) ^ 
a(pp) + a{pn) ' 

where a = d 2 a/dMdy\ y=0 , as a function of y/r at Pbcam = 450 GeV/c, corresponding 
to the proposed CERN experiment 0.Q The parton distributions are the same as in 
Fig. 3. Because u > d in the proton, the asymmetry is positive for sets with d — u zero 
or small at large x. The asymmetry is reduced and can even become negative for sets 
with large d — u, in particular the CTEQ1M set. A measurement of this asymmetry 
in the \fr = 0.1 — 0.3 range to an accuracy of a few percent will therefore provide 
a powerful discriminator of the possible behaviours of the d — u difference. Note 
also that the sensitivity of the asymmetry (12) to d — u increases at forward rapidity 
(equivalently xp > 0): at large y, where uu annihilation again dominates, we have 
simply that A DY ~ (u — d)/(u + d), which is equivalent to (|TT| ) with 1 — 2Z/A = 1. 

In conclusion, we find that deep inelastic data do not unambiguously determine 
u and d. However, if a reasonable behaviour is imposed on d — u at small x (mo- 
tivated by a small breaking of meson exchange degeneracy) then the data indicate 

2 The predictions for the corresponding Fermilab energy are virtually indistinguishable. 
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a small flavour asymmetry with d > u for x < 0.01. Global fits to data imply that 
the d > u inequality persists to larger x, but it is possible to introduce more flexible 
parametrizations of d — u. We show that a comparison of Drell-Yan production on 
hydrogen and deuterium targets will provide a good determination of d — u in the 
medium x range. Finally, it is worth noting that a small, positive d — u is a common 



feature of most model calculations JIT], |T2| of the sea quark flavour asymmetry, in 



particular those involving mesonic contributions to the deep inelastic structure func- 



tions. For example, the Drell-Yan predictions of the chiral quark model of ref. |12 
are very similar to those of our standard Dq fit. 
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Figure Captions 



[1] The parton distributions of the proton at Q 2 = 20 GeV 2 corresponding to (a) 
set D' of MRS @ and to (b) set CTEQ1M of g . 

[2] The upper plot shows the accumulated contribution to the Gottfried sum rule, 
(6), as a function of x, the lower limit of integration, which is obtained from 
sets of partons with differing behaviour of d — u, as shown in the lower plot (for 
Q 2 = 7 GeV 2 ). The shaded bands contain a family of parton solutions with 
d — u given by (8) with —8 < 7 < 32 and 5 = 0; the central curves in the bands 
corresponds to set D of @] with 7 = 0. The dashed curves correspond to a fit 
with 5 7^ in (8), and the dot-dashed curves correspond to set CTEQ1M of [Q. 
The central part of the figure shows the integrand of the Gottfried sum rule, 
together with NMC data || which we have corrected, at small x, for deuteron 
screening effects. 

[3] The ratio of dilepton yields per nucleon from tungsten and isoscalar (deuterium 
and carbon) targets as a function of a; (target). The data are from ref. J7|; 
the open circles at small x show the ratio before the correction for nuclear 
shadowing. The curves are the full next-to-leading order QCD predictions using 
the same parton sets as in Fig. 2. Thus the shaded band corresponds to a range 
of distributions with d—u given by (8) with —8 < 7 < 32 and 5 = 0; the central 
curve in the band corresponds to set Dq of with 7 = 0. The dashed curve 
corresponds to a fit with 5 = 300, 7 = —60 in (8), while the dot-dashed curve 
corresponds to set CTEQ1M of 0. 

[4] Predictions for the Drell-Yan asymmetry Ady = (cr(pp) — cr(pn))/(a(pp) + 
a{pn)), where a = d 2 a/dMdy\ y=0 , as a function of y/r = M/y/s with pbeam = 
450 GeV/c. The curves are based on next-to-leading order calculations using 
the same parton distributions as in Fig. 3. 
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